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On the dynamics of large- O (A^)-symmetric quantum systems at finite temperature 
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Time evolution of a perturbed thermal state is studied in a quantum-mechanical system with 
O (N) symmetry. In the limit of large A^, time dependence of O (A'')-singlet expectation values can 
be described by classical equations of motion in a one-dimensional potential well. Time depen- 
dence of the perturbation is then described by a linear differential equation with time-dependent 
periodic coefficient. This equation, depending on the parameters, admits either exponentially grow- 
ing/decaying or periodically oscillating solutions. It is demonstrated that only the latter possibility 
is actually realized, thus in such a system there is no redistribution of initial perturbation over all 
A'^ degrees of freedom. 
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This work was motivated by several closely related 
nonstationary problems in quantum field theory. The 
first two are Schwinger pair creation in a strong electric 
field [1] and decay of cosmological constant in de Sitter 
space [2, 3, 4]. Another problem is "fast scrambling" of 
information in matrix quantum mechanics considered re- 
cently in [5, ()]. In all cases initial state of the system 
is classically stable but decays due to quantum effects. 
Typically, only a few degrees of freedom out of some large 
or infinite number are excited initially, but due to inter- 
actions the energy and information content of the initial 
state is redistributed more or less equally among all de- 
grees of freedom. 

In the context of such problems, one is usu- 
ally interested in expectation values of the form 
(m| e~'^*Oe'^* |m), where H is the Hamiltonian and 
O is some quantity which characterizes a perturbation 
of the system (for example, field strength for the decay 
of strong electric field). Practical calculations arc usu- 
ally done either using the perturbation theory, semiclas- 
sical approximation or numerical methods. It is clear 
that perturbative expansion can only lead to the re- 
sults of the form (m| e-'(^"+3-^^)*C'e'(^«+*'-^0* |m) = 
{in\ O \in) + aigt + a2g^t^ + which are valid either for 
small t ov g. It is not surprising that such a dependence 
on time was indeed obtained, for example, in perturba- 
tive calculations of graviton loops in de Sitter background 
[3] . On the other hand, semiclassical calculations similar 
to that of Schwinger [1] usually fail to take into account 
back-reaction processes. 

The aim of this Letter is to obtain an exact descrip- 
tion of the time evolution of a perturbed thermal state 
in a model system with very large number N of interact- 
ing degrees of freedom. The Hamiltonian of this model 
system is O (A'^)-symmetric and has the following form 



[7, 8, 9]: 
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where (j)a, tTo are the coordinate and momentum opera- 
tors and V (x) grows monotonously for a; > and does 
not depend on N . The terms in V {x) which are pro- 
portional to correspond to interactions involving k 
degrees of freedom. Note that each degree of freedom 
interacts with each other degree of freedom with equal 
strength. In the limit of infinite N such system can be 
thought of as a sort of mean-field approximation to a 
quantum theory of scalar field. Indeed, if one consid- 
ers scalar fields on A^-dimensional lattice, each lattice 
site interacts equally with N neighbors, and the limit 
N 00 leads to the standard mean-field approxima- 
tion. This analogy closely follows the discussion of "fast 
scrambling" in [5] , where the following problem was con- 
sidered: suppose one has a system with a large number 

of degrees of freedom. How should these degrees of 
freedom interact in order to redistribute perturbations of 
a small number of degrees of freedom over all the system 
most efficiently and quickly? The answer suggested in [■^)] 
is precisely that each degree of freedom should interact 
with each other degree of freedom with equal strength. 
Thus one might hope that the system described by the 
Hamiltonian (1) is a "fast scrambler". 

Here the following particular problem will be con- 
sidered: suppose the system (1) is initially in a state 
of thermal equilibrium with density matrix po ~ 

cxp (^—(3H^, where (3 is the inverse temperature. 
Then the system is perturbed in such a way that the 
density matrix pq ((j)a, ipa) the coordinate representa- 
tion is shifted by a constant vector S^a- Pin{4>a,4>'a) = 

N 

Pa i4>a ~ £,a,(t>'a- ia)- '^^ assumcd to be finite in 

a=l 

the limit N 00. After that the system is allowed to 
evolve freely. What will be the time dependence of the 

expectation value (j)a (t) = Tr I pin<j>a (t) ), where (j>a (t) 
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is the coordinate operator in the Heisenberg representa- 
tion? Will it decay with time because the energy of an 
initial perturbation is redistributed over all A'^ degrees 
of freedom, or will it oscillate periodically? The large- 

analysis presented here suggests that only the latter 
possibility is actually realized for any potential V (x). 

Heisenberg equations of motion for (pa and tTq are 

^(^a = TTa, ^ TTa = -2V' (x) ^a, where X (t) = 

jf ^ (j^a (t) and V (x) is the derivative of the potential 

a=l 

over X. Differentiating once again over time, one obtains: 



N 



N 

Y,^l~W'{x)x = 



N 



H ~W {x) ~ MV' {x) 



(2) 



(3) 



where in the last line the operator ^ tt^ was expressed 

in terms of the Hamiltonian (1). 

A crucial step now is to take the expectation values of 
the r.h.s. and l.h.s. of (2) and (3) over the initial state 
with the density matrix ■ For the ground state of the 
Hamiltonian (1), or, more generally, for the thermal state 
po in the limit of large A^, expectation values of multiple 
powers of a; factorize [7]: lim Tr (poi*^) = Tr {p^xf" . 

Since the density matrix pm was obtained from po by a 
shift of the coordinates 4>a, it is easy to show that the 



lim Tr (p,-. 



same property holds also for p^, 

Tr {pinx)^ . Factorization property is usually considered 
only for ground states of large- TV theories [7, 9], however, 
it is easy to show that it holds also for the initial thermal 
state Po and for the perturbed state at t > 0. To this end 

one can consider the path integral for Tr {j)in4>a [t]^ = 

Tr {4>a (0) cxp {i J dr (^H + Ja (r) 0a) ) ) and note that 

it has a saddle point which dominates in the large-A^ 
limit. The existence of such a saddle point in the path in- 
tegral corresponds to factorization property in the canon- 
ical formalism. Here the integral goes over the T-shaped 
contour in the complex plane which consists of the dou- 
ble line of length t on the real axis and of the line of 
length (3 on the imaginary axis. At the ends of this 
line the usual periodic boundary conditions are imposed. 
At the intersection of the two lines, a singular source 
Ja (t) = —ia ^ 6 (r, 0) is inserted into the path integral 
in order to reproduce the initial perturbation. Path in- 
tegral over the imaginary axis yields then the thermal 
ground state and path integral over real time describes 
evolution of the perturbed state in time. However, the 
equations of motion (2) and (3) arc most easily analyzed 
in the canonical formalism. 

Assuming factorization, introducing the expectation 

values 4>a {t) = Tr (^Pin(i>a , x{t) = Tr (pi„x (t)), and 
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FIG. 1: A typical shape of the potential V (x). 



taking into account that the Hamiltonian (1) does not 
depend on time, one can obtain a closed system of equa- 
tions for (pa (t), X (t): 



d^<Pa 
dt^ 



-2c^aV' [X) 

d'^x d ~ 



V{x)=4x [Vix)--Ti- {p,nH 



Initial conditions for the equations (4) and (5) are: 



<l^a{0)=^a, ^(0)=0, 

^(0) = xo + |, §(0) = 



(4) 
(5) 
(6) 



(7) 



where = Tr {pox) is the expectation value of x 
in the thermal state po. The last initial condition 
can be justified in the following way. For any quan- 

tum state l^*) one has (^l ^ |*) = J2 (*l ^aT^a + 

a=l 

N / _ _ X 

^a<Pa I*) = i E / d^'^a i^^a^^^ + (0a*) = 

a— 1 ^ ^ 

N , _ _ X 

^ E jd^cpa V^<l^a^^^ ~ {cpa^) gf^* • Thc CXpCCta- 
a— 1 ^ ^ 

tion value (^| ^ x j^*) is equal to the sum of the above 
expression over all eigenstates \'^n) of the Hamiltonian 
(I). The states with the wave functions and are 
both the eigenstates of the Hamiltonian (1) with the same 
energy, and therefore both of them enter the thermal den- 
sity matrix with equal weights. This means that for each 
pair \E'„ and 4'„ the terms in the brackets give zero in 
total. Therefore ^ (0) = Tr (p„^ = 0. 

Consider now the potential V {x). It depends on 
the initial conditions due to the term Tr I pmH ) = 



Tr ypin^j + NV{x{0)), which can be expressed in 
terms of xq and ^a- To this end note that Tr (^Pm-f-^ — 
Tr (^Po-f'^ and make use of the virial theorem, which 
states that Tr (po^) = Tr (pgi^^j = NxqV (xq). 
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Thus the potential V (x) is: 

V (x) = AxV (x) - Ax {xoV (xo) + V{x (0))) (8) 

It is easy to show that if V (x) is a monotonic con- 
tinuous function for a; > 0, equation (5) always de- 
scribes a finite periodic motion in a one-dimensional po- 
tential well. Namely, there is always such Xf < x(0), 
X/ > that V{xf) = V{x{Q)) a.ndV{x) < V {x (O)) for 
Xf < X < X (0), so that X (t) moves periodically between 
Xf and X (0) with some period T. A typical shape ofV (x) 
is illustrated on Fig. 1. Note that V (x) is not bounded 
from below, but depends on x (0) in such a way that the 
motion of x is always finite and is always restricted to 
X > 0. 

The equations (4), (5) together with the initial condi- 
tions (7) are the main result of this work. Thus in the 
large A'' limit x moves classically in the one-dimensional 
potential well V {x) , while (t) is a solution of the lin- 
ear differential equation with time-dependent coefhcicnt 
2V' (x). Since x (t) is periodic, one can readily apply the 
Floquet theorem, which states that the equation (4) with 
periodic x (t) can have either two independent oscillating 
solutions or two solutions one of which grows and the 
other decays exponentially. In the first case the eigenval- 
ues of the Floquet matrix have unit absolute values and 
are complex conjugate to each other: |Ai| = IA2I = 1, 
Ai = A2. In the other case, which corresponds to para- 
metric resonance between x and (f>a, the eigenvalues Ai, 
A2 are real and A1A2 = 1. 

One might therefore ask, whether by tuning the pa- 
rameters of the system it is possible to achieve such a 
parametric resonance between x (t) and 4>a (t)? The an- 
swer to this question seems to be negative. Exponentially 
growing (j)a (t) is clearly prohibited on general physical 
grounds, since in this case any small perturbation causes 
the energy of all N degrees of freedom to be collected into 
one mode. While such solutions may be in principle rel- 
evant for spontaneous symmetry breaking, they are very 
unlikely for monotonic potentials V (x). I was not able to 
find a general proof of this fact, however, in numerical in- 
vestigations of the equations (4), (5) for several different 
V (x) no regions of parametric resonance were found. 

But probably there can be isolated points in parame- 
ter space, for which (j>a (t) decays exponentially? At first 
sight on equations (5), (4) and (7), it seems that such 
possibility may indeed be realized, and the system (1) 
may indeed "scramble" an initial perturbation. However, 
this turns out to be impossible as well. The reason for 
that is, at the bottom of the fact, the invariance of the 
system (f ) under time reflection. Indeed, (j)a (t) can only 
decay exponentially if the initial state with (j)a (0) ~ ^a, 
jf. ipa (0) = is an eigenstate of the Floquet matrix of 
the equation (4) with eigenvalue A < I. Otherwise (pa (t) 
will also contain the mixture of the second independent 
solution of (4), which grows exponentially. Therefore the 
time derivative of (j)a should be equal to zero at times 
t = nT: ^ (j)a (nT) = 0. Since x (t) undergoes one- 




FIG. 2: The solution (j>a (t) of (4) (projected on the direction 
of <j)a (0)) for the potential V (x) = + ^ a.t the inverse 
temperature /3 = 0.5. Dots denote the numerical result for 
0a (t) obtained by solving the saddle-point equations for (jta (t) 
and X (t) in the path integral representation. 



dimensional motion in a potential well V (x) with the 
initial condition ^ x (0) = 0, x (t) is symmetric under 
time reversal: x (t) ~ x {nT — t) for any n. The equa- 
tion (4) is therefore also invariant under time reversal, 
and if one evolves the system forward in time from t = 
to t ^ nT and back to t = 0, one should arrive at the 
initial value 4>a (0). On the other hand, this value should 
be equal to X^^ipa (0). We have thus to conclude that if 
A is real, it can only be equal to one, and the equation 
(5) has no exponentially decaying solutions. There can 
only be the solutions which oscillate periodically with 
time. For example, a typical solution of (4) for the po- 

tcntial Vix) = ^ + ^,Nv(f)=^ + ^ 

at the inverse temperature (3 = 0.5 is plotted on Fig. 2. 
Dots denote numerical solution for (jja (t), which was ob- 
tained by numerically solving the saddle-point equations 
for X (t) and (pa {t) [7, S, 9] on the T-shaped contour in the 
plane of complex time, as explained above. At not very 
large t the amplitude of oscillations of (j)a (t) indeed de- 
creases, but at larger times 0a (t) oscillates in a beat-like 
manner. Note that such solution is completely different 
from the classical one, where both x (t) and (t) oscil- 
late with equal frequencies without any beats. The beats 
in (j)a (t) are therefore a purely quantum phenomenon. 

It appears that the system described by the Hamilto- 
nian (1) does not have an asymptotic property of redis- 
tributing a perturbation of a small number of degrees of 
freedom over all configuration space at any finite temper- 
ature. Instead there is a continuous swapping of energy 
between an excited degree of freedom and all other de- 
grees of freedom. This may be related to the fact that 
the system (I) has as many integrals of motion (the com- 
ponents of angular momentum) as there are degrees of 
freedom. Therefore the interactions between all A^ de- 
grees of freedom should be arranged in some more com- 
plicated way in order to achieve "scrambling" . Perhaps 
one should indeed give them a structure of large matrices, 
as suggested in [-5, G]. However, the equations (4) and (5), 
which describe the dynamics of O (A'^)-symmetric large- 
N model and not just the properties of its ground state. 
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may be of interest by themselves. It might be interesting 
to find similar equations in higher-dimensional large- 
theories or in matrix models. 
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